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A substantial fraction of the energy density of the universe 
may consist of quintessence in the form of a slowly- rolling 
scalar field. Since the energy density of the scalar field gen- 
erally decreases more slowly than the matter energy density, 
it appears that the ratio of the two densities must be set to 
a special, infinitesimal value in the early universe in order 
to have the two densities nearly coincide today. Recently, we 
introduced the notion of tracker fields to avoid this initial con- 
ditions problem. In the paper, we address the following ques- 
tions: What is the general condition to have tracker fields? 
What is the relation between the matter energy density and 
the equation-of-state of the universe imposed by tracker so- 
lutions? And, can tracker solutions explain why quintessence 
is becoming important today rather than during the early 
universe? 



I. INTRODUCTION 

Quintessence^ has been proposed as the missing en- 
ergy component that must be added to the baryonic an 
matter density in order to reach the critical densityB 
Quintessence is a dynamical, slowly-evolving, spatially 
inhomogeneous component with negative pressure. An 
example i^a. scalar field Q slowly rolling down its poten- 
tial V(Q)BoE!l For a general scalar field, the pressure is 
Pq = \Q 2 — V, and the energy density is pq — \Q 2 + V. 
For a slowly rolling scalar field, the pressure can be neg- 
ative if the kinetic energy is less than the potential en- 
ergy. For quintessence, the equation-of-state, defined as 
wq = Pq/pq, lies between and — 1. Depending on 
V(Q), wq may bejconstant, slowly- varying, rapidly vary- 
ing or oscillatory.El 

A key problem with the quintessence proposal is ex- 
plaining why pq and the matter energy density should be 
comparable today. There are two aspects to this prob- 
lem. First of all, throughout the history of the universe, 
the two densities decrease at different rates, so it ap- 
pears that the conditions in the early universe have to 
be set very carefully in order for the energy densities to 
be comparable today. We refer to this—issue of initial 
conditions as the "coincidence problem. "E3 This is a gen- 
eralization of the_flatness problem described by Dicke and 
Peebles in 1979;li3 the very same issue arises with a cos- 
mological constant as well. A second aspect, which we 



call the "fine-tuning problem," is that the value of the 
quintessence energy density (or vacuum energy or cur- 
vature) is very tiny compared to typical particle physics 
scales. The fine-tuning condition is forced by direct mea- 
surements; however, the initial conditions or coincidence 
problem depends on the theoretical candidate for the 
missing energy. 

Recently, we introduced a form of quintessence called 
"tracker fields" which avoids the coincidence problcm.Ell 
Tracker fields have an equation-of-motion with attractor- 
like solutions in which a very wide range of initial condi- 
tions rapidly converge to a common, cosmic evolutionary 
track. The initial value of pq can vary by nearly 100 
orders of magnitude without altering the cosmic history. 
The acceptable initial conditions include the natural pos- 
sibility of equipartition after inflation - nearly equal en- 
ergy density in Q as in the other 100-1000 degrees of 
freedom (e.g., Qq, ps 10~ 3 ). Furthermore, the result- 
ing cosmology has desirable properties. The equation-of- 
state wq varies according to the background equation- 
of-state wb- When the universe is radiation-dominated 
(wb = 1/3), then wq is less than or equal to 1/3 and pq 
decreases less rapidly than the radiation density. When 
the universe is matter-dominated (wg — 0), then wq is 
less than zero and pq decreases less rapidly than the mat- 
ter density. Eventually, pq surpasses the matter density 
and becomes the dominant component. At this point, 
Q slows to a crawl and wq — > — 1 as J1q — > 1 and the 
universe is driven into an accelerating phase. These prop- 
erties seem to match current observations well.E£l 

In some cases, wq is very nearly equal to wb during 
the radiation- and matter-dominated epochs, and some 
authors have focused on the fact that the Q energy den- 
sity decreases ai—nearly the same rate as the background 
energy density.tL3'tl3 However, our focus is on the fact 
that the cosmology is independent of initial conditions, 
which includes but does not require wq ~ wb during the 
matter- and radiation-dominated epochs. Hence, our use 
of the term "tracker" is meant to refer to solutions join- 
ing a common evolutionary track, as opposed to tracking 
closely the background energy density. 

It is also an interesting point that some tracker solu- 
tions do not require small niaas. parameters to obtain a 
small energy density today;oliH whether this is a satis- 
factory solution to the fine tuning problem is debatable, 
though. We do not explore this issue in this paper. 

The tracker solution differs from a classical dynam- 
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ics attractor because it is time-dependent: during track- 
ing, Qq increases steadily. This contrasts with the solu- 
tions based on exponential potentials recently discussed 
by Ferreira and Joycet2l in which Hq is constant during 
the matter dominated epoch. For this kind of solution, 
Qq is constrained to be small < 0.15) at the be- 
ginning of matter domination in order that large-scale 
structure be formed, but then it cannot change there- 
after. Hence, it remains a small, subdominant compo- 
nent. (See, however, discussion of undershoot and over- 
shoot solutions later in this paper for a possible way 
around this constraint.) The tracker solution is more de- 
sirable, though, because it enables the Q _ener gy to even- 
tually overtake the matter density and induce a period of 
accelerated expansion, which produces a cosmology more 
consistent with measurements of the matter density, large 
scale structure, and supernovae observations. 1513 

An important consequence of the tracker solutions is 
the prediction of a relation between wq and £Iq today.tfl 
Because tracker solutions are insensitive to initial condi- 
tions, both wq and Qq only depend on V(Q). Hence, for 
any given V(Q), once S1q is measured, wq is determined. 
In general, the closer that £Iq is to unity, the closer wq 
is to —1. However, since Vt m > 0.2 today, there is a suf- 
ficient gap between Q,q and unity that wq cannot be so 
close to —1. We find that wq > —0.8 for practical mod- 
els. This wq-^Iq relation, which makes the tracker field 
proposal distinguishable from the cosmological constant, 
will be explored further in this paper. 

The purpose of the present paper is to expand on our 
introductory article on tracker fields and the coincidence 
problem. In particular, we want to go beyond the specific 
examples of tracker potentials studied before and address 
in a general way the following questions (the relevant 
section is shown in parentheses): 

• What is a tracker solution? (II. B) 

• What potentials V(Q) have tracker solutions and 
which do not? (III.A-C,F) 

• How does convergence of diverse initial conditions 
to a common track occur? (III.D) 

• What additional conditions on tracking potentials 
are required to make them useful and practical? 
(III.H) 

• How does tracking result in the prediction of an 
Q,q-wq relation? (IV) 

• What range of wq today is possible for tracking 
solutions according to the wq-^Iq relation? (IV) 

• Why is the Q-field first beginning to dominate at 
this late stage of the universe rather than at some 
early stage? (V) 



II. DEFINITIONS AND BASIC EQUATIONS 

We shall consider a scalar field with present equation- 
of-state — 1 < wq < in a flat cosmological background 
(consistent with inflation) . The ratio of the energy den- 
sity to the critical density today is Qq for the Q-field 
and fi TO for the baryonic and dark matter density where 
Q m + Qq = 1. We use dimensionless units where the 
Planck mass is M v = 1. 



A. The Tracker Equation 

The cquation-of-motion for the Q-ficld is 

Q + 3HQ + V = (1) 

where 

H 2 = Q = n{ Pm + Pr + ig 2 + V) (2) 

where a is the Robertson- Walker scale factor, p m is the 
matter density, p r is the radiation energy density, and 
K = The definition of the equation-of-state is 
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It is extremely useful to combine these relations into 
an unfamiliar form, which is the form we would like the 
reader to have in mind when we refer to the "equation- 
of-motion" : 
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where x = (1 + wq)/(1 — wq) = \Q 2 jV is the ratio of 
the kinetic to potential energy density for Q and prime 
means derivative with respect to Q. The ± sign depends 
on whether V' > or V' < 0, respectively. The tracking 
solution (to which general solutions converge) has the 
property that wq is nearly constant and lies between wb 
and -1. For 1 + w Q = 0(1), Q 2 w Q Q H 2 and the 
equation-of- motion, Eq. (||), dictates that 
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for a tracking solution; we shall refer to this as the 
"tracker condition." 

An important function is Y = V"V/ (V) 2 , whose prop- 
erties determine whether tracking solutions exist. Taking 
the derivative of the equation-of-motion with respect to 
Q and combining with the equation-of-motion itself, we 
obtain the equation: 
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where x = d Inx/d In a and 'x = d 2 \a.x/d In a 2 . We will 
refer to this equation as the "tracker equation." 



B. What is a tracker solution? 

We shall use the term "converging" to refer to 
the attractor-like behavior in which solutions to the 
equation-of-motion are drawn towards a common solu- 
tion, Q(t). We shall argue that converging solutions with 
wq < wb solve the coincidence problem. A tracker solu- 
tion is one which undergoes long periods of converging. 
The converging period should be long enough so that any 
residue of initial condition is washed away and the solu- 
tion today only depends on parameters in the potential. 
For a potential V(Q) — M v(Q/M) (where v is a dimcn- 
sionless function of Q/M), there is a family of tracker 
solutions parameterized by M. The value of M is deter- 
mined by the measured value of £l m today (assuming a 
flat universe). 

We make the (fine) distinction between "tracking" 
and "converging" because tracker solutions typically go 
through some periods when solutions are not approach- 
ing one another. In particular, during transitions in the 
background equation-of-state, such as the transition from 
radiation- to matter-domination, there are brief intervals 
during which wq changes rapidly and general solutions 
are not drawn towards the tracker solution; fortunately, 
these transitory periods are too short to spoil the advan- 
tages of tracking solutions. 

To be complete in classifying the possibilities, one 
should also consider potentials in which Q first passes 
through a regime where solutions converge and then ends 
up in a regime where they do not. We refer to these as 
"hybrid models" below. By the time the second regime is 
reached, sufficient convergence of solutions may have al- 
ready occurred to insure that the desired range of initial 
conditions track along a common solution in the second 
regime. 



III. PROPERTIES OF TRACKING SOLUTIONS 
AND TRACKING POTENTIALS 

In the following section, we prove the key properties of 
tracker solutions. Each subsection leads with a summary 
statement of the result for those who wish to skip the 
detailed mathematical arguments. 



A. What potentials have converging behavior and 
tracker solutions and which do not? 



Our central theorem is that tracking behavior with 
wq < wb occurs for any potential in which T = 
V"V/(V') 2 > 1 and is nearly constant (\d(T-l)/Hdt\ < 
|T — 1|) over the range of plausible initial Q. For the 
cases of interest, the range extends to V(Q) equal to 
the initial background energy density pb down to V(Q) 
equal to the background density at matter-radiation 
equality, a span of over 100 orders of magnitude. Al- 
though wq < wb is what is desired for most practi- 
cal applications, we also show that tracking behavior 
with wb < wq < (1/2)(1 + wb) is possible provided 
1 — (1 — wb)/{6 + 2wb) < T < 1 and nearly constant. 
Tracking does not occur for T < 1 — (1 — wb)/(6 + 2u>b). 

Hence, testing for the existence of tracking solutions 
reduces to a simple condition on V(Q) without having to 
solve the equation-of-motion directly. In particular, the 
condition that T w constant can be evaluated either by 
testing 
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over range of Q corresponding to the allowed initial con- 
ditions; the middle expression is easily computed know- 
ing V(Q) only without having to solve an equation- 
of-motion. (Here we have used the tracker condition 
Q/H w v/Qq sw V'/V which applies for the tracker so- 
lution.) An equivalent condition is that that Ar/r <C 1, 
where Ar is the difference between the maximum and 
minimum values of T over the same range of Q. The con- 
dition r > 1 is equivalent to the constraint that \V /V\ 
be decreasing as V decreases. These conditions encom- 
pass an extremely broad range of potentials, including 
inverse power-law potentials (V(Q) = M 4+a /Q a for a > 
0) and combinations of inverse power-law terms (e.g., 
V(Q) = M 4 exp(M/Q)). Some potentials of this form 
are suggested by particle physics models with d^namipal 
symmetry breaking or nonperturbative effects O'EHcil 

If wq < wb, converging behavior does not occur for 
potentials in which strictly increases as V de- 

creases (r < 1). This category includes quadratic po- 
tentials and most examples of quintessence models in the 
literature. Instead of converging behavior, these models 
require specially tuned initial conditions to obtain an ac- 
ceptable value of Oq today, as discussed in Section III.B. 

For (1/2) (1 + wb) > wq > wb, converging behavior 
does occur if V'/V is strictly increasing as V decreases 
(r < 1); however, as we explain in Section III.F, poten- 
tials with wq > wb do not produce viable cosmologi- 
cal models: if wq > wb and Qq > 1/2 today (as sug- 
gested by current observation) then J7q must exceed O m 
in the past and there is no period of matter-domination 
or structure formation. 
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Hybrid potentials are possible in which converging be- 
havior only occurs for a finite period during the early uni- 
verse, provided that the time is sufficiently long to bring 
together solutions whose initial conditions span the range 
of practical interest. For example, one can construct po- 
tentials in which decreases at first as Q rolls down- 
hill, and then |V'/V| begins to increase. Another possi- 
bility is that |F'/F| increases at first and wq > wb, as 
discussed in the previous paragraph. We pointed out that 
this condition allows convergence, but cannot be main- 
tained up to the present for practical reasons of structure 
formation. However, the condition could be maintained 
for a long, finite period in the early universe so long as it 
is terminated before matter-domination. 



B. Why do models with increasing | V'/ | and 
wq < wb fail to solve the coincidence problem? 

For potentials in which | V'/V| increases as V decreases 
(r < 1), the LHS of the the equation-of- motion Eq. (|J) is 
increasing. If wq < wb, 1/ \J£Iq on the RHS is decreas- 
ing as Q rolls downhill. Hence, the tracker condition, 
|V'/V| ~ 1/ •v/Oq, cannot be maintained and wq cannot 
be maintained at a nearly constant value (different from 
— 1) for any extended period. 

In particular, extrapolating backwards in time, | V'/V| 
is strictly decreasing and 1/ WCIq is strictly increasing. 
The only way to satisfy the equation-of- motion, Eq. (Ji|), 
is to have wq approach — 1. But, this corresponds to 
the Q-field freezing at some value Q = Qi after only a 
few Hubble times. The energy density of the frozen field 
is only slightly higher than the current energy density. 
Consequently, to obtain this cosmological solution, one 
has to have special initial conditions in the early uni- 
verse (after inflation, say) that set Q precisely to Qi and 
the Q-energy density to pQ i , a value nearly 100 orders 
of magnitude smaller than the background energy den- 
sity. It is precisely this tuning of initial conditions (the 
coincidence problem) which we seek to avoid. 

We note that many quintessence models (wq < 0) 
discussed in the literature fall into this non-tracking 
class and require extraordinary tuning of initial condi- 
tions. Simple examples include the harmonic potential, 
V(Q) = M 2 Q 2 , and the sinusoidal potential, V(Q) = 
M 4 [cos (Q/ f) + 1]. In these models, the initial value of 
Q must be set to a specific value initially in order to 
obtain the measured value of £Iq today. 



C. What determines the tracker solution? 

Our claim is that solutions converge to a tracker solu- 
tion if: (a) r(Q) = V"V/(V') 2 is nearly constant; and (b) 
T > 1 for w Q < w B or T < 1 for (1/2)(1 + w B ) > w Q > 
wb- We have explained in the previous section why T 



must be greater than unity if wq <wgor else the track- 
ing condition, V'/V ~ 1/ WSIq, cannot be maintained. A 
similar argument requires 1 — (1 — wb)/(6+2wb) < T < 1 
if (1/2) (1 + wb) > wq > wb- But, as we shall see, we 
also need that T be nearly constant in order to have con- 
verging behavior. 

If r is nearly constant, then the tracker equation, 
Eq. (^|), implies that there is a solution in which wq is 
nearly constant and x and its time-derivatives are negli- 
gible. In this case, we have that the equation-of-state for 
the Q-deld is nearly constant: 



wq 



w B - 2(r - i) 
i + 2(r-i) ' 



(8) 



If r > 1, the value of wq must be less than % to sat- 
isfy the tracker equation, which means that the Q-energy 
red shifts more slowly than the background energy. For 
inverse power-law potentials V(Q) — M A+a /Q a , T = 
1 + a -1 , and the relation .ate have derived matches the 
relation in our first paperJi-1 For V(Q) = M 4 exp(l/Q), 
r = 1 + 2Q, which satisfies our condition for Q <C 1. 
For r < 1, the same relation says that the value of 
wq must be greater than wb] but we also require that 
r > 1 — (1 — wb)/(6 + 2wb) in order to satisfy the crite- 
rion, wq < (1/2)(1 + wb), necessary to have converging 
solutions (see Section III.F). 

If r varies significantly with Q, then one can find a 
wide and continuous distribution of Q, J7q, and wq which 
satisfies the tracker equation. Hence, there is no single 
solution to which solutions converge. 

In sum, we have succeeded in expressing our conditions 
for tracking solutions in terms of T, which depends en- 
tirely on the functional form of V. A simple computation 
of V"V / (V 1 ) 2 determines if V admits tracking solutions 
or not. The case of interest for quintessence is wq < wb, 
in which case the condition for a tracker solution is T > 1 
and nearly constant. 



D. How is the tracker solution approached beginning 
from different initial conditions? 



We shall call the tracker solution Q(t) and the energy 
density of the tracker solution as a function of time Pq (t) . 
This subsection explains in rough detail how solutions 
converge to the tracker solution for any initial pq that 
lies between the initial background energy density, psi, 
and the current critical density, p C Q. This range stretches 
nearly 100 orders of magnitude. In particular, we explain 
how the convergence to the tracker solution is different 
if initially p B > Pq > Pq versus p C Q < p Q < pq. For 
simplicity, we confine ourselves to the case wq < wb ■ We 
complete the proof of convergence in the next subsection. 

The equation-of- motion, Eq. (0) can be rearranged as: 
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3 v / k(1 + w Q ) 



l = A(i)-l; (9) 



where we have restricted ourselves for simplicity to po- 
tentials with V < 0. Then, since V < 0, the RHS of 
this equation is a balance between a positive semi-definite 
term (A) and a negative term. The tracker condition 
corresponds to the two conditions: 1 — Wq significantly 
different from zero and wq nearly constant. The latter 
requires near balance on the RHS above (A ss 1) so that 
d hix/d In a is nearly zero. 
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FIG. 1. Energy density versus red shift for the evolution 
of a tracker field. For computational convenience, z — 10 12 
(rather than inflation) has been arbitrarily chosen as the ini- 
tial time. The white bar on left represents the range of initial 
Pq what leads to undershoot and the grey bar represents over- 
shoot, combining for a span of more than 100 orders of magni- 
tude if we extrapolate back to inflation. The solid black circle 
represents the unique initial condition required if the missing 
energy is vacuum energy density. The solid thick curve repre- 
sents an "overshoot" in which pq begins from a value greater 
than the tracker solution value, decreases rapidly and freezes, 
and eventually joins the tracker solution. 

First, consider the "overshoot" case in which pQ is 
initially much greater than the tracker value pQ. For 
simplicity, let us assume that Q is released from rest. 
The evolution goes through four stages, illustrated in 
Figures 1-3. The potential for this exampleis V(Q) — 
M /Q 6 . (For this and all subsequent figures, the choice 
of z = 10 12 has been chosen as the initial time for con- 
venience of computation and illustration; a realistic fig- 
ure would have initial z corresponding to the inflationary 
scale.) 

1. V /V and S!q are so big initially that A 3> 1. So, 
1 d In x 2 d In wq 



+ 1. This means that Q becomes large and V de- 
creases very rapidly as Q runs downhill. 

As Q runs downhill, V /V and are decreasing. 
Consequently, A begins to decrease and ultimately 
reaches a value of order unity, one of the require- 
ments for a tracker solution. However, wq has been 
driven towards +1; up to this point, the RHS of 
Eq. (^) has been positive, so there has been no op- 
portunity for wq to decrease. As a result, 1 — Wq is 
too small, or, more specifically, the kinetic energy is 
too large for Q to join the tracker solution. Hence, 
Q rolls farther down the potential, overshooting the 
tracker solution. 

Once the tracker solution is overshot, A becomes 
less than unity and the RHS of Eq. @ changes sign. 
wq now decreases from +1 towards —1. One might 
wonder what happens when wq crosses through the 
tracker value; why doesn't Q track at this point? 
The answer is that there is now the problem that 
A is too small. So, the RHS of Eq. (§) remains too 
negative and wq continues to decreases and heads 
towards —1. 

Once wq reaches close to —1, Q is essentially frozen 
at some value Qf and, consequently, V and V /V 
are frozen. However, A is now increasing since Qq 
is increasing - even though pq is nearly constant, 
Pb is decreasing. As A increases to order unity, the 
sign of Eq. (||) changes once again; wq increases 
from — 1; the field runs downhill; and, the sign of 
Eq. (^) changes yet again. After a few oscillations, 
the terms in Eq. (|^) settle into near balance and Q 
is on track. 
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and 



wq 



is driven towards its maximal value, 



wq 
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FIG. 2. A plot of wq vs. red shift for the overshoot solu- 
tion shown in Figure 1. wq rushes immediately towards +1 
and Q becomes kinetic energy dominated. The field freezes 
and wq rushes towards —1. Finally, when Q rejoins the 
tracker solution, wq increases, briefly oscillates and settles 
into the tracker value. 
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FIG. 3. A plot of i/6 = (1/6) d Inx/d In a for the over- 
shoot solution shown in Figures 1 and 2. At late times when 
Q settles into the tracker solution, x is small and wq is nearly 
constant. During the overshoot phase, x undergoes large pos- 
itive and negative changes, as described in the text. 

Next, consider the "undershoot" case in which pq is 
initially much less than the tracker value pQ and Q is 

released from rest. This corresponds to Q ^> Q initially. 
By assumption, V and | are much smaller than 

the tracker value. Consequently, 1/ \JQ,q is larger than 
the tracker value. The only way to satisfy the equation- 
of- motion, Eq. (Q) is for wq to approach —1 so that the 
coefficient of 1/ ^/Hq is nearly zero. This condition corre- 
sponds to a very small kinetic energy density or Q nearly 
constant. Hence, the field remains nearly "frozen", and 
V and V'/V are nearly constant as the universe evolves. 
The situation is identically to beginning with Step 4 
above, and the scenario proceeds just as described there. 

In sum, the field either drops precipitously past the 
tracker value and is frozen (overshoot), or it begins with 
a value less than the tracker solution (undershoot) and is 
frozen. In either case, it proceeds from the frozen state 
to joining the tracker solution. In the case of undershoot, 
the frozen value Qf is simply the initial value of Q. For 
the overshoot case, Q begins by going through a kinetic 
energy dominated period in which Q cx a~ 3 . If the initial 
PQi < PBi, where is the background radiation density, 
then a cx i 1 / 2 and the frozen value of Q is 

Qf~Qi + JJjkh (11) 



where the subscript refers to the initial values of Q and 
£Iq. (If initially pQi > psi, then a cx t 1 / 3 , and 

Q/«Q 4 + /f (1 + Jln^); (12) 

V 47T 2 p B 

in this case, Qf is so large that Q remains frozen up to the 
present time.) For the overshoot case, Qi is typically very 
small compared to unity and Qqi — 0(1). Consequently, 
the frozen value Q / depends on QQi only. 

Initial conditions in which Qi is non-zero do not change 
the story significantly. If Qi is very large, then the ini- 
tial behavior is kinetic energy dominated, and the evolu- 
tion proceeds similar to the overshoot case. Initial fluc- 
tuations in Q also do not change the story since they 
are exponentially suppressed once the potential becomes 
non-negligible and the field is driven towards the tracker 
solution (see Appendix). 

The possibility of overshoot and undershoot allows a 
new possibility for the case of exponential potentials re- 
cently discussed by Ferreira and Joyce.EEl The exponen- 
tial potential is a special example of a tracker solution 
in which J7q is constant during the matter dominated 
epoch. The practical problem with this model, as noted 
by Ferreira and Joyce, is Qq is constrained to be small 
(TIq < 0.15). At the beginning of matter domination, 
Oq must be small in order that large-scale structure be 
formed; but then it cannot change thereafter. Hence, it 
remains a small, subdominant component. This argu- 
ment presumes, however, that Q is already on track at 
the beginning of matter domination. It is possible to tune 
initial conditions so that Q overshoots the tracker solu- 
tion initially and does not join the tracker solution until 
just very recently (red shift z = 1). Then, the constraint 
on VLq is lifted. 

E. What are the constraints on the initial value of Q 
and pQ? 

Suppose that the tracker solution corresponds to Q = 
Q today and the Q has converged to a tracker solution. 
Then, whether pq is initially smaller than the tracker 
value and frozen at some Q = Qf equal to its initial value, 
or pq is initially larger than the tracker value and falls 
to Qf « ^SftQi/AiT, it is necessary that Qf be less than 
Qo in order that the field be tracking today. This is not a 
very strong constraint. Since Qo = 0(1) for most track- 
ing potentials, this only requires that p c o < PQi < pBi 
initially. Hence, initial conditions in which Q dominates 
the radiation and matter density are disallowed because 
Q falls so fast and drops to such a low point on the poten- 
tial that it has not yet begun tracking today. However, 
initial conditions in which there is rough equipartition 
between pq and the background energy density are al- 
lowed, as well as initial values of pq ranging as low as 
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100 orders of magnitude smaller, comparable to the cur- 
rent matter density. The allowed range is impressive and 
spans the most physically likely possibilities. 



F. Is the tracker solution stable? 

What has been shown so far is that, whether the ini- 
tial conditions correspond to undershoot or overshoot, 
Q soon reaches some frozen value Qf which depends on 
the initial Qi. Then, after some evolution, Qq increases 
to the point where |V"/V| ~ 1/tJQq and, according to 
the equation-of-motion, wq moves away from — 1 and the 
field begins to roll. What remains to be shown is that 
solutions with wq not equal to the tracker solution value 
converge to the tracker solution. Or, cquivalently, we 
need to show that the tracker solution is stable. 

Now, consider a solution in which wq differs from the 
tracker solution value wq by an amount 8. Then, the 
master equation can be expanded to lowest order in 5 
and its derivatives to obtain after some algebra: 

5 + 3[^Ob + 1) - w Q }5 + |(1 + w B )(l - w Q )S = (13) 

where the dot means d/d In a as in the tracker equation. 
The solution of this equation is 



(14) 



where 

3 1 

1 = -^(wb + 1)-w ] 



±% 2 V 18(1 + WB ^ ~ ™ o) ~ 9 ^ WB + w °l 2 - ( 15 ) 

The real part of the exponent 7 is negative for wq be- 
tween — 1 and wo = \ ( 1 + wb ) , which includes our entire 
range of interest. So, without imposing any further con- 
ditions, this means that 8 decays exponentially and the 
solution approaches the tracker solution. As 8 decays, it 
also oscillates with a frequency described by the second 
term. See Figure 4. 
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FIG. 4. The convergence of different initial conditions 
to the tracker solution. As derived in the text, wq decays 
exponentially fast to the tracker solution combined with small 
oscillations. All the curves are for V(Q) = Af 4 /Q 6 . The solid 
curve is the overshoot case from Figure 1. The thin dashed 
curve with w ~ is the tracker solution which is overlaid for 
most z by the dash-dotted curve, which represents a slightly 
undershooting solution. 

In deriving Eq. ([To]), we have assumed that T is 
strictly constant, independent of Q, which is exactly 
true for pure inverse power-law (V ~ 1/Q a ) or expo- 
nential (V ~ exp(/3Q)) potentials. The same result 
holds if \d(T - \)/Hdt\/\T - 1| < 1 {i.e., T - 1 varies 
with Q but only by a modest amount) over the plausi- 
ble range of initial conditions ranging from V(Q) f=a pbi 
to V(Q) w p eq (where pBi is the initial background en- 
ergy density after inflation, say, and p eq is the energy 
density at matter-radiation equality). The condition is 
equivalent to |r'/[r (V'/V)} \ < 1. In this limit, T and 
the tracker value of wq change adiabatically as Q rolls 
downhill, satisfying the tracker equation with x being 
negligibly small, as discussed in the Appendix. The con- 
straints on wq are the same as above. However, some im- 
portant differences from the constant T case are pointed 
out in the last section. 

Throughout most of our discussion in this paper, we 
have considered the case wq < wb- However, our con- 
vergence criterion, wo < |(1 + wg), includes wq > wb 
or, equivalently, 1 — (1 — wb)/(6 + Ppb) < T < 1, 
as also found by Liddlc and Scherrer.Ej An example is 
V ~ Q a with a > (6 + 2w B )/(l ~ w B ); a > 6 for 
wb = and a > 10 for wb = 1/3. Let us suppose 
we reached the present Qq after tracking down this po- 
tential. Because wq > wb in these potentials, it must 
be that f2g exceeds FIb extrapolating backwards in time. 
Consequently, there is no period of matter-domination 
or structure formation, and these models have no practi- 
cal interest. However, see the discussion of hybrid models 
below for a variation on these models that may be viable. 
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G. Borderline models and hybrid models 

For completeness, we consider two special classes of 
potentials, borderline trackers in which T = 1 and hybrid 
models in which T > 1 at first and then r < 1 . 

The borderline case corresponds to V ec-oxiiWQ), 
which has been studied by several authorsBEltiJllil For 
this case, the tracker equation for x = x = Q demands 
that wq — wb and, therefore, JIq is constant. Hence, for 
this case, the tracker solution corresponds to maintain- 
ing a constant ratio of quintessence to background energy 
density. The only deviation occurs during the transition 
from radiation- to matter-domination when x becomes 
non-negligible, but this is a small effect. 

Because fig is constant throughout the matter- 
dominated epoch, these models have limited practical 
utility. £Iq must be small (< 15%) at the onset of matter- 
domination in order not to disrupt structure formation. 
(Quintessence suppresses the growth rate.) But, then, 
since Qq is constant, SIq remains small forever. Conse- 
quently, the models require il m > 85*2^ inconsistent with 
a number of determinations of mass,Ej and the universe 
never enters a period of accelerated expansion, incon- 
sistent with recent measurements of thfi luminosity-red 
shift relation for Type IA supernovae.Ej (The overshoot 
scenario may lift the Q m > 85% constraint but, as dis- 
cussed in Section D, introduces fine tuning which defeats 
the whole purpose of the scenario.) 

Hybrid models have the property that solutions con- 
verge to a tracker solution at the early phase of evolu- 
tion but cease to converge after a certain point due to a 
change in the shape of the potential as Q rolls downhill. 
One can imagine a sufficiently long convergence regime 
that all or most plausible initial conditions have collapsed 
to a common tracker solution before the second regime 
begins. Effectively, this has the desired feature that a 
wide range of initial conditions lead to the same final 
condition. The models may be somewhat artificial in 
that the current cosmology is very sensitive to where the 
transition occurs. For example, consider the case where 
wq < wb but r undergoes a transition from T > 1 (con- 
verging) to T < 1. Recall that T < 1 corresponds to 
|^'/^l increasing as V decreases. We have shown that, 
extrapolating backwards only a small interval in time, 
the field Q must have been frozen at a value not so dif- 
ferent from the current value (assuming the field is rolling 
today and wq < 0). In this kind of hybrid model, the 
transition to V < 1 must be set so that the transition 
occurs so that Q is near the frozen value, which requires 
delicate tuning of parameters in the potential. 

A different example is where 1 — (1 — wb)/{6 + 2wb) < 
r < 1 and wq > wb during the early stages of the 
universe. We have argued that these conditions pro- 
duce converging behavior but, if the conditions continue 
to the present, there is no period of matter-domination 



or structure formation (see Section III.F). However, one 
can imagine hybrid models in which these conditions 
[1 — (1 — wb)/{6 + 2wb)] < T < 1 and wq > wb are 
satisfied for some period early in the history of the uni- 
verse, providing a finite period of converging behavior. 
Later, as Q moves down the potential, V changes form so 
that wq < wb ■ Viable cosmological models of this type 
can be constructed in which Qq does not dominate the 
universe during the matter-dominated epoch until near 
the present time. 

H. Some additional practical considerations 

We have discovered a wide class of potentials that ex- 
hibit tracking behavior. This guarantees that a wide 
range of initial conditions converge to a common tracking 
solution, but the convergence may take longer than the 
age of the universe in some cases. In particular, if one 
assume equipartition after inflation, say, and the initial 
V is too far above the tracker solution, then Q falls pre- 
cipitously, overshoots the tracker solution, and freezes at 
some Q = Qf. For some potentials satisfying the condi- 
tion r w constant > 1, Qf may be so large that the field 
does not begin to roll and track by the present epoch. If 
Q just started to roll by the present epoch, then it would 
behave exactly as a cosmological constant until now, and 
so the model is trivially equivalent to a A model. As a 
practical consideration, we demand that a field starting 
from equipartition initial conditions should start rolling 
by matter-domination, say, so that the model is non- 
trivial. This imposes a mild added constraint on poten- 
tials, V(Q). For this purpose, rough estimates suffice. 

Equipartition at the end of inflation, when there are 
hundreds or perhaps thousands of degrees of freedom 
in the cosmological fluid, means that the Q-field has 
Qi 10~ 3 . Beginning from equipartition, Q falls to some 
value Qf where it freezes. According to the equation-of- 
motion, Q remains frozen until •v/Oq cx \/V / H increases 
to where V'/V ~ I/^/Qq or, equivalently, 

Yl {Qi)>H \ Zeq ) = ^-^^. (16) 

For V(Q) cx 1/Q a , this imposes the constraint 

a 2 1 
Qf 2 -n Q°al q > (17) 

where Qo is the present value of Q. Since today 

v£ H ° -*th' (18) 

we obtain Qq ~ ^figa; 2 . From Eq. (pi]), we also know 
that 
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47T o 

— q) ~ ~ iQ- 



(19) 



Combining the above relations one gets the restriction on 
a 



a 2 > 



io 12 (j^) (1+ */ 2 



(20) 



where we have taken a e? ~ 10~ 4 . This approximate rela- 
tion leads to a > 5. Figure 5 confirms this result showing 
that the a = 1 model starts rolling much later than equal- 
ity. So, if one restricts to pure inverse power-law (\/Q a ) 
potentials, our constraint that Q begin from equiparti- 
tion and roll before matter-radiation equality constrains 
us to a > 5. 




Z+1 

FIG. 5. A comparison of the overshoot for three different 
models beginning from f2j = 10 -3 (equipartition). The thick 
solid line is for V(Q) = M 4 /Q 6 , the thick dash-dotted line 
is for V = M 4 [exp(l/Q) — 1], and the thick long dashed line 
is for V(Q) = M 4 /Q. In all three examples, Q falls rapidly 
downhill and freezes. In the first and second examples, Q 
begins to roll again and joins the tracker solution before mat- 
ter-radiation equality; the third example, which violates the 
condition derived in the text, does not begin to roll again by 
the present epoch. 



IV. THE Q,q-Wq RELATION 

An extremely important aspect of tracker solutions is 
the Qq-wq relation, or, equivalcntly, the £l m -WQ relation 
which it forces. For any given V — M A v(Q /M) (where 
v is a dimensionless function of Q/M), Q and Q are to- 
tally determined independent of initial conditions by the 
tracker solution. The only degree of freedom is the M- 
parameter in the potential. M can be fixed by imposing 
the constraint that the universe is flat and Q m = 1 — £Iq is 
determined by measurement. There is, then, no freedom 
left to independently vary wq. This is the explanation 
of the Qq-wq relation, a new prediction that arises from 



tracker fields. The Qq-wq relation is not unique because 
there remains the freedom to change the functional form 
of V(Q). Even so, wq is sufficiently constrained as to be 
cosmologically interesting. 

The general trend is that wq -1 as £Iq — ► 1. The 
fact that fl m > 0.2 observationally means that f2g < 0.8 
and wq cannot be very close to —1. How small wq can 
be is model-dependent. We are most interested in the 
smallest values of wq possible since the difference from 
— 1 determines how difficult it is to distinguish the track- 
ing field candidate for missing energy from cosmological 
constant. 

In Figure 6, we show the £Iq-wq relation for a series of 
pure, inverse power-law potentials, V(Q) oc 1/Q a . The 
general trend is that wq increases as a decreases. The 
constraint given at the end of the previous section is that 
a > 5 (in order that Q be rolling by matter-radiation 
equality beginning from equipartition initial conditions) . 
For J7q = 0.8, the smallest value of wq is -0.52, which 
occurs for a = 5. This is a large difference from — 1 
obtained for a cosmological constant. 

However, a lower value of wq can be easily achieved for 
a more generic potential with a mixture of inverse power- 
laws, e.g., V{Q) = M 4 exp(l/Q). For these models, Q 
is small initially. If the potential is expanded in inverse 
powers of Q, 1/Q a , then it is dominated in the early 
stages by the high-a terms. Hence, the effective value of 
a is much greater than 5 before matter-radiation equal- 
ity, and we easily satisfy the constraint that the field be 
rolling before matter-radiation equality beginning from 
equipartition initial conditions. On the other hand, the 
value of Q at the present epoch is large, and the potential 
is dominated by the a ~ 1 terms in its expansion. Con- 
sequently, wq can be even lower today than in the pure 
power-law case. For J7q = 0.8, we obtain wq — —0.72 
for the exponential potential, which is in better-accord 
with recent constraints on wq from supernovae.u 

It is difficult to go below this limit without artificially 
tuning potentials unless we relax our constraints. For 
example, consider the highly contrived potential V(Q) = 
+ (Jiff) hi which we have intentionally chosen ex- 
ponents differing by six orders of magnitude in order to 
obtain a small wq today. The second term in the poten- 
tial dominates before equality and insures that the field 
is rolling by that point in time. The first term dominates 
at late times and makes the equation of state wq very low 
(because this term is very flat). For this example, we find 
wq = —0.98 for S!q = 0.8. However, we had to choose 
a pair of terms with exponents differing by six orders of 
magnitude. The lesson of this exercise (and related tests) 
is that the exponential potential is a reliable estimate for 
the minimal wq possible for generic, untuned potentials. 
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FIG. 6. The £Iq-wq relation for various potentials assum- 
ing a flat universe tt m = 1 — Qq, where wq represents the 
present value of wq. The potentials and notation are the 
same as in Figure 5. 

In Figure 6, we illustrate the Uq-wq relation for sev- 
eral potentials. Only the V ~ 1/Q 6 and V ~ exp(l/Q) 
potentials satisfy the condition that Q is rolling by 
matter-radiation equality beginning from equipartition 
initial conditions. In Figure 7, we illustrate the effec- 
tive Wq that would be measured using supernovae or 
cosmic microwave background measurements, using the 
V ~ exp(l/Q) potential as defining the boundary of 
minimal values of wq possible for the tracker field case. 
This boundary assumes that we satisfy the strict condi- 
tion that the field be rolling by matter-radiation equal- 
ity beginning from equipartition initial conditions. If 
we relax this condition and allow a somewhat narrower 
range of initial conditions, then general potentials of the 
form V ~ T,Ck/Q k (such as V ~ 1/Q) are allowed and 
wq can be somewhat smaller (see Figure 6). Hence, in 
Figure 7, the boundary can relax somewhat downward 
(dashed line) but it is difficult to obtain wq < —0.8 or 
Wq < —0.75. Because wq is evolving at recent times, 
the value obtained from measurements at moderate to 
deep red shift will differ from the current value shown in 
Figure 6. For tracker potentials, the effect of integrat- 
ing back in time over varying wq turns out to be well- 
mimicked by a model with constant wq — Wq^ that has 
the same conformal distance to last scattering surface. 
For the case J7 m = 0.2, for example, Figure 6 shows that 
wq = —0.72 today, but Figure 7 shows that the measured 
wq would be w e ff — —0.63. 

The result is exciting because the VLq-wq relation 
and the constraint on TO > 0.2 today creates a suffi- 
ciently large gap between Wq* and —1 that the tracker 
candidate for missing energy should be distinguishable 
from the cosmological constant in near-future cosmic mi- 
crowave background and supernovae measurements. 
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FIG. 7. A plot of Wq ! versus £l m = 1 — Qq, showing 
the minimum Wq' possible for tracker solutions. The solid 
line is the lower boundary assuming the constraint that the 
Q-field begins with equipartition initial conditions and be- 
gins rolling before matter-radiation equality. The dashed line 
is the lower boundary if this condition is relaxed to allow gen- 
eral V ~ T,Ch/Q k . As explained in the text, Wq ' is the value 
that would be measured in supernovae and microwave back- 
ground experiments which effectively integrate over a varying 
Wq. 



V. WHY IS THE UNIVERSE ACCELERATING 
TODAY? 

We have proven in this paper that tracker potentials 
resolve the coincidence problem for quintessence. For a 
very wide range of initial conditions, cosmic evolution 
converges to a common track. The tracker models are 
similar to inflation in that they funnel a diverse range of 
initial conditions into a common final state. The models 
have only one important free parameter (M) which is 
fixed by the measured Qq = 1 — Q m . 

Some of the mathematical properties of tracking so- 
lutions have been noted before. £ox_fixp(/3Q) and pure, 
inverse power-law potentialsoQ'oEil However, the ex- 
traordinary insensitivity to initial conditions and the po- 
tential application to the coincidence problem was not 
explored. The present work is important because it 
shows that the properties are shared by a much wider 
class of more generic potentials. "Generic potentials" in- 
clude, for example, all V"s which can be expanded as 
a finite or infinite sum of terms with inverse powers of 
Q, which is much more general than the special cases 
of a single inverse-power or a pure exponential. We use 
V = Af 4 exp(l/Q) as an example of this more generic 
class, although our conclusions would remain the same 
for more general V — T,Ck/Q k - 

Extending the tracker behavior to generic potentials 
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may be important because, as we shall argue below, 
they have properties not shared by the special cases 
(V ~ 1/Q a and V ~ exp(/?Q)) that make plausible why 
Qq only begins to to dominate and initiate a period of 
accelerated expansion late in the history of the universe. 

To understand this conclusion, we need to change our 
approach. Up to this point in the paper, we have imag- 
ined fixing M to guarantee that Qq = 1 — f2 m has the 
measured value today. This amounts to considering one 
tracker solution for each V(Q). Now we want to consider 
the entire family of tracker solutions for each given V(Q) 
and consider whether Qq is more likely to dominate late 
in the universe for one V(Q) or another. 
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FIG. 8. A plot of P/Pq versus t, where Q,q oc t p and 
Po is the initial value of P. The plot compares pure inverse 
power-law (V ~ 1/Q a ) potentials for which P is constant 
with a generic potential (e.g., V ~ exp(l/Q)) for which P 
increases with time. 



In general, Qq is proportional to a 3 ^ WB WQ ' oc 

t 2(w B -w Q )/(l+w B )^ where W(J have shown in Eq. (|) that 



W B - WQ 



2(T-l)(tfl B + l) 



i + 2(r-i) 

Hence, we find Qq oc t p where 



P 



4(r-i) 
i + 2(r-i)' 



(21) 



(22) 



A satisfying explanation as to why Qq dominates at 
late times, rather than early, would be if Qq automati- 
cally changes behavior at late times. This is not the case 
for the two special cases (V ~ 1/Q Q and V ~ exp(j3Q)). 
As the universe evolves, T — 1 and, hence, P are con- 
stant; consequently, Qq grows as the same function of 
time throughout the radiation- and matter-dominated 
epochs. However, these are the exception, rather than 
the rule. 



For more general potentials, P increases as the uni- 
verse ages. Consider first a potential which is the sum of 
two inverse power- law terms with exponents a% < a.2- 
The term with the larger power is dominant at early 
times when Q is small, but the term with the smaller 
power dominates at late times as Q rolls downhill and 
obtains a larger value. Hence, the effective value of a 
decreases and r — 1 oc 1/a increases; the result is that P 
increases at late times. For more general potentials, such 
as V ~ exp(l/<5), the effective value of a decreases con- 
tinuously and P increases with time. Figure 8 illustrates 
the comparison in the growth of P. 

How does this explain why Qq dominates late in the 
universe? Because an increasing P means that Qq grows 
more rapidly as the universe ages. Figure 9 compares 
a tracker solution for a pure inverse power-law poten- 
tial (V ~ 1/Q 6 ) model with a tracker solution for V ~ 
exp(l/Q), where the two solutions have been chosen to 
begin at the same value of Qq . (The start time has been 
chosen arbitrarily at z = 10 17 for the purposes of this 
illustration.) Following each curve to the right, there is a 
dramatic (10 orders of magnitude) difference between the 
time when the first solution (solid line) meets the back- 
ground density versus the second solution (dot-dashed 
line). That is, beginning from the same Qq, the first 
tracker solution dominates well before matter-radiation 
equality and the second (generic) example dominates well 
after matter-domination. The difference is less dramatic 
as a increases for the pure inverse power-law model and 
becomes negligible for a > 15. Of course, the model 
appears more contrived. But, more importantly, as a 
increases, the value of wq today (given Q m > 0.2) ap- 
proaches zero and the universe does not enter a period 
of acceleration by the present epoch. Hence, a signifi- 
cant conclusion is that the pure exponential and inverse 
power-law models are atypical; the generic potential has 
properties that make it much more plausible that Qq 
dominates late in the history of the universe and induces 
a recent period of accelerated expansion. 
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FIG. 9. A plot comparing two tracker solutions for the 
case of a V ~ 1/Q 6 potential (solid line) and al/~ exp(l/Q) 
potential (dot-dash). The dashed line is the background den- 
sity. The two tracker solutions were chosen to have the same 
energy density initially. The tracker solution for the generic 
example (V ~ exp(l/Q)) reaches the background density 
much later than for the pure inverse-power law potential. 
Hence, Qq is more likely to dominate late in the history of 
the universe in the generic case. 

In sum, the general tracker behavior shown in this pa- 
per goes a long way towards resolving two key issues: 
the coincidence (or initial conditions) problem and why 
Pq is dominating today rather than at some early epoch. 
And, it leads to a new prediction - a relation between 
Q m = 1 — Qq and wq today that makes tracker fields 
distinguishable from a cosmological constant. 
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VI. APPENDIX 

In this Appendix, we discuss the convergence to the 
tracker solution when Y = V"V/ (V) 2 varies with Q. We 
wish to show that convergence occurs if the variation of 
r over the plausible range of initial conditions (varying of 
100 orders of magnitude in energy density) is nearly con- 
stant. An example is V ~ exp(l/Q) for which T = 1 + 2Q 
and Q <C 1 for the plausible range of initial conditions. 

The condition that Y is nearly constant means pre- 
cisely that d(Y — 1)/Hdt <C (Y — 1) or, equivalently, 
|r'/[r(V/V)]| < 1. In this case, Y is nearly constant 
over a Hubble time. Hence, we can consider an adiabatic 
approximation for the tracker solution in which Qq and 
wq satisfy the tracker equation Eq. (^) with x and x negli- 
gibly small. Suppose wq and Q are both perturbed from 



this tracker solution by amounts Sw and SQ(t). From the 
definition of wq, we have that 

5w = (l-w )—SQ-(l + w )—SQ. (23) 
Po Po 

From the cquation-of-motion Eq. (Q) we know that 
Po 



,-3(wq-w b ) 



^0 ^Tq 

Consequenty, Eq. (|3|) implies 

Sw ~ exp [— 3/2 (ws — wq) In a] SQ 



(24) 



(25) 



In particular, this equation shows that Sw decays if SQ 
decays. 

To show the SQ decays, we start from the more stan- 
dard form of the equation-of-motion: 



Q + 3HQ + V' = 0, 
and obtain the perturbed equation 

SQ + 3HSQ + VqSQ = 0. 
Changing the variable to dr = Hdt, we obtain 



d 2 SQ _ d fdSQ 
dr 2 ~ ~Hdt VWdt 



1 » dSQ 



(26) 
(27) 

(28) 



where m = (Ht)- 1 = 3{w B + l)/2. Eq. @ then be- 
comes 



d 2 SQ 
dr 2 



+ ( 3- m) fU^0. 



where 



C 



T r/l 

H 2 



r(! - w Q )(2 + w B +w q ) 



(29) 



(30) 



in the adiabatic approximation. The solution to Eq. ( |29| ) 

is SQ — A exp (/3r) , where 



(3 2 + (3 - m)(3 + C = 



with solutions 
where 



D 
~2 



1 i 1/1 



4C 
B 2 



B = 3 — m = —(1 — %). 



(31) 



(32) 



(33) 



Note that B > and C > 0. Hence, we have that SQ ex- 
ponentially decays, and, by the argument that preceded, 
Sw cx SQ also decays exponentially. Combining our re- 
lations for SQ and Sw, we can reproduce the result in 
Eq. (pa) obtained for the constant Y case. 
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If Q has spatial fluctuations, Eq. ( p9| ) must be modified 
by a positive term proportional to k 2 5Q. The effect is to 
increase C and modify the oscillation frequency. How- 
ever, the exponential suppression of the fluctuations is 
retained once the field starts approaching the attractor 
solution. Hence, even if the initial conditions result in 
significant fluctuations after the field is frozen (in either 
the undershoot or overshoot case) , the initial fluctuations 
are erased as Q converges to the tracker solution. 
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